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We introduce a simple model to describe the frictional properties of granular media under shear. 
We model the friction force in terms of the horizontal velocity x and the vertical position z of the 
slider, interpreting z as a constitutive variable characterizing the contact. Dilatancy is shown to 
play an essential role in the dynamics, inducing a stick-slip instability at low velocity. We compute 
the phase diagram, analyze numerically the model for a wide range of parameters and compare our 
results with experiments on dry and wet granular media, obtaining a good agreement. In particular, 
we reproduce the hysteretic velocity dependence of the frictional force. 



I. INTRODUCTION 

Problems related to interfacial friction are very important from a practical and conceptual point of view , and in 
spite of its wide domain of application sliding friction is still not completely understood. In order to construct efficient 
machines in engineering science, or to understand geophysical events like earthquakes, it is necessary to understand 
several aspects of friction dynamics. Beside usual solid-solid contacts, the sliding interface can be lubricated with 
molecular fluids or filled by a granular gauge, and the problem becomes strongly dependent on the internal dynamics 
of the material itself. In experiments on lubricated surfaces with thin layers of molecular fluids the friction force 
depends on the thermodynamic phase of the lubricant, which depends on its turn on the shear stress [Q. The friction 
force is thus directly related to the microscopic dynamics of the system and a description of sliding friction cannot 
be achieved without a good microscopic understanding of the problem. Understanding the frictional properties of 
granular matter turns out to be an even harder task, since basic problems like stress propagation in a static packing 
remain largely unsolved due to the disordered nature of the stress repartition inside the medium. Moreover, when a 
granular medium is sheared, it reorganizes modifing the geometrical disorder. The microscopic arrangement of the 
grains and their compaction have an important effect on the friction since, in order to deform the medium one has to 
overcome several geometrical constraints. 

The understanding of sheared granular media has recently advanced thanks to experiments [^-^ and numerical 
studies §-§. The response to an external shear stress can be characterized by the dilatancy, which measures the 
modification of the compaction of the granular medium during the flow 0. We note that in both lubricated and 
granular interfaces the friction force has a dynamical origin. Since a sheared material modifies its own internal state 
fiuidizing or changing structure, a natural approach to the problem is to describe phenomenologically this change of 
state and to relate it to the macroscopic friction force. As we discussed previously, a complete theoretical description 
of sheared granular media is still not available, so that the analysis should strongly rely on experimental data. 

Recent experiments focusing on the stick-slip instability induced by friction in sheared granular layers, helped to 
elucidate the role of compaction and the microscopic origin of slip events. In particular, accurate measurement of the 
friction force and of the horizontal and vertical positions of the slider Q have permitted to emphasize the connections 
between dilatancy and friction. The apparatus used was composed by a slider dragged at constant velocity by a spring 
whose elongation measured the applied shear stress. The surface of the slider was roughened in order to avoid slip at 
the surface of the medium and so that friction would crucially depend on the internal structure of the medium. At 
low velocity, a stick slip instability was observed and related to the modification in the granular compaction. 

Friction of granular layers has been mainly studied in the framework of geophysical research p|,pO[| using rate and 
state constitutive equations [pdl list where the friction force is a function of an auxiliary variable describing the state 
of the interface. In this approach, one assumes that the microscopic events causing the movement of the slider are 
self-averaging and neglects the fiuctuations. The quantities used in the constitutive equations are thus mean field-like. 
This assumption should be valid for sliding friction experiments on granular materials, where the size of the grain is 
much smaller than the length of the slider, so that the variables used in the model (velocity, displacement or friction 
force), are well defined macroscopical quantities. 

The constitutive variable, related to the microscopic dynamics of the system, describes the dynamical history of 
the interface. In the case of solid-solid interfaces this variable was associated with the age of the contacts and 
described two opposite effects: the age of the contact increases the static friction force and the displacement of the 
slider renews the interface continuously so that the friction force decreases with velocity. Lubricated systems have 
been approached similarly using the rate of fiuidization as a constitutive variable which captures two different 

effects. On one hand the confinement of the thin fiuid layer induces a glassy transition resulting in a large static friction 
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force. On the other hand an appHed shear stress increases the temperature of the medium, favoring fluidization, thus 
reducing the friction force, which crucially depends on the ratio between the strength of the two effects. 

In the case of granular media, a parameter suitable to characterize the frictional behavior is the compaction of 
the layers or the height of the slider which can be measured experimentally. Also in this case we can identify the 
competition between two opposite effects: the velocity of the slider keeps the layer dilated, lowering the friction force, 
and the weight of the slider induces recompaction. In this paper we present a model which includes these two effects 
in the framework of rate and state constitutive equations to describe typical effects like the stick-slip instability or 
the force- velocity hysteretic loop. 

In Sec. H we concentrate on the description of the model, in Sec. [II we describe the main results obtained b 
numerical integration of the model, in Sec. IV we present a stability analysis and the phase diagram. Finally, Sec. 
presents a discussion and a summary of our results. 



II. THE MODEL 



Here we write rate and state constitutive equations in order to describe the frictional properties of granular media. 
The dynamics of the sliding plate is described by two constitutive equations. The first one is simply the equation 
of motion for the slider block driven by a spring of stiffness k and submitted to a frictional force, which depends 
on velocity and dilatancy. The second equation is the evolution law for an auxiliary variable characterizing the 
dilatancy, which we will identify with the vertical position of the slider. This model could in principle be applied to 
geophysical situations, although in that case instead of a single elastic constant k, strain is mediated via the material 
bulk elasticity. 

The frictional properties of a granular medium depend explicitly on its density: a dense granular medium submitted 
to a tangential stress tends to dilate, i.e to modify the granular packing and thus the friction force. It is not simple 
to measure granular density especially for non homogeneous systems, but global changes can be characterized by 
the vertical position of the sliding plate, which is thus an excellent candidate to describe the state of the system. 
Therefore, in agreement with Ref. [Q, we write the equation of motion for the slider block as 

mit ^ k{Vt - x) ~ F{z,x), (1) 

where m is the mass of the sliding plate, x its position, k the spring constant, V the drag velocity, and F{z,x) the 
friction force depending on the velocity x and on the height of the plate z. 

If the slider is at rest, we need to apply a minimal constant force Fq in order for it to move. When the force 
exceeds Fq, the slider moves and dilation will occur, reducing the friction. When the layer is fully dilated the friction 
force reduces to Fd < Fq. We assume that the friction force is velocity dependent when the layer is partially dilated 
(z < Zjn), and becomes independent on velocity in the stationary state, when the granular medium is fully dilated 

— ^m) • 

In summary (in the case i > 0), we write the friction force as 

F(z, x)=Fa- (3^-^ - (2) 

The first two terms in Eq. (2) give the friction force at rest {x = 0) as function of z. In the fully expanded phase 
{z = Zjn), the friction term is F = < Fq, while in the compacted phase Fq = Fd + (3zm/R {z — 0). The velocity 
dependence is linear, mediated by the factor z — z,„ which vanishes when the bed is fully dilated. These equations 
should be compared with those presented in Ref. 0|, where the second term in Eq. (2) is not present. 

In Eq. (2) F(z, x) depends explicitly on 2, which describes the vertical displacement of the slider. In order to 
complete the description of the dynamics, we must specify the evolution equation for z. We write the law controlling 
the dilation of the granular medium during shear as 

In Eq. (3) the second term dilates the support and can be seen as the response of the granular medium to the 
external tangential stress: when submitted to a shear rate i, the medium dilates and z increases. The factor (z — z^) 
reduces to zero when the bed is fully dilated and Zm can be identified with the maximal height. 

The first term allows for recompaction under the slider weight: in the case a; = the plate falls exponentially fast. 
At high velocity this term will not perturb significantly the system and the dynamics will be stationary. We are 
interested in the small velocity limit: Eqs. (2,3) will display an instability below a critical drag velocity Vc, as we will 
show in Sec. 
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It is useful to rewrite the system of equations in terms of dimensionless variables 

^^K' ^^"^v' ^"i?' ^"""r' 

m^m-.^ V^V—, v = v—, F, = —, (3=—. 

Defining I — Vt — x — Fd, the system becomes 

i^V-d, (4) 
rhv = I + {z - z~rri){v + (3), (5) 

z = -- - (z - z™){;. (6) 
Assuming that these equation are valid for i > 0, we can analyze them for different spring constants, velocity. 

III. NUMERICAL SIMULATIONS 

We numerically solve the model (Eqs. (4-6)) using the fourth order Runge-Kutta method and assuming that the 
slider plate sticks when its velocity is zero. We concentrate our analysis on two sets of parameters. The first set 
corresponds to experiments carried out with a dry granular medium. We compute the typical force-velocity diagram 
in order to fix the parameters. Then using the same parameters we test the validity of our model calculating other 
quantities such as the slider velocity during a slip event, the spring elongation or the vertical displacement. 

A second set of parameters is used to model wet granular media. We recover the instability at low velocity and low 
spring force and study the evolution of dilatancy and spring elongation before reaching the steady-state. 

A. Dry granular media 

Dry granular media exhibit stick-slip instabilities for relatively high velocities and it is difhcult to achieve complete 
vertical displacement of the slider. For this reason the steady sliding regime has not been studied in detail in 
experiments. In order to quantitatively test our model we adjust the parameters to fit the experimental results. We 
present in Fig. |l| the force- velocity curve during slip comparing the experimental data from Ref. ||^ with the result of 
the integration of the model. The parameters used are given in the caption. The model is able to accurately describe 
the first part of the hysteretic loop (when the velocity increases) , but slight deviations appear for small velocities for 
which also the experimental uncertainties are larger. 

We numerically integrate the model using the previously obtained parameters, varying the spring constant and 
the driving velocity. For slow velocity and a small spring constant the system exhibits typical stick-slip dynamics. 
Fig. ^ shows the evolution of the variables of the model in this case: the first plot (Fig. 2(a)) shows the variation 
of the spring length which decreases abruptly at a regular frequency, when the horizontal plate position increases 
(Fig. 2(b)). Fig. 2(c) represents the velocity of the plate which is followed by an increase of the vertical position of 
the plate (Fig. 2(d)). We show in Fig. || a more detailed study of the slider velocity during a slip event. Near the 
transition between the stick-slip and the steady sliding regime the slider velocity appears to be almost independent 
on the driving velocity, in agreement with experiments. The stick-slip instability of the model is ruled by Eqs. (4-6) 



and we present in Sec. ^ the dynamical phase diagram computed by a linear stability analysis. When the slider is 
driven slowly the energy injected into the granular medium cannot keep the layers dilated and the motion stops after 
a short change in the horizontal position (slip event). 

If we increase V ot k, the energy induced by the shear is sufficient to maintain the granular layer dilated and 
the system evolves to a steady sliding state (cf. Fig. Q), which is stable with respect to small perturbations. This 



stationary state, corresponds to a stable fixed point of Eqs. (4-6) (see Sec. iV). If the drag velocity is very large 
the steady sliding state becomes unstable due to inertial effects (m 7^ 0) and the slider oscillates harmonically with 
frequency uj = ^/kjm. This effect was experimentally observed in Ref. ||] . We have plotted the result in Fig. || for two 
different perturbations, in order to show that the amplitude of the cycle depends on the strength of the perturbation. 
A typical measurement performed in the framework of geophysical research 10|, is the variation of the friction 



force with respect to a rapid change of the driving velocity. We have simulated this effect, and we show the result 
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m Fig. |. An increase of the driving velocity is followed by an increase of the friction force which then relaxes to a 
smaller value. 

The phase diagram corresponding to the three different dynamical behaviors can be calculated analytically. We 
present the result in Sec. [V , where we study the linear stability of the model. 



B. Wet granular media 



The analysis performed in Sec. Ill A can be repeated for wet granular materials. The dynamics in this case is more 
stable and the stick-slip regime is more difhcult to obtain experimentally, since the instability occurs at very slow 
velocity. In the wet case, the presence of water changes the dynamics of the grains. Under shear, grains reorganize 
submitted to the fluid viscosity but here we neglect the small hydrodynamic effects and consider only the grain 
dynamics with suitable parameters. Using the new class of parameters, we solve numerically Eqs. (4-6) and identify 
two regimes: steady sliding at high V and stick-slip instability otherwise (see Sec. IV for more details). In Fig. |^ we 
show a typical plot of the different quantities in the stick-slip regime. The period of the oscillations is bigger than 
in the dry case, and the fluctuations of the elongation smaller. One of the main difference with the dry case is the 
value of 7] which governs the relaxation process and which is greater in the wet case as an effect of immersion. In 
Fig. 1^ we show the steady state found at high velocity. It is interesting to remark that this behavior can be perfectly 
recovered with a simplified model, presented in Ref. which however does not give rise to stick-slip instabilities. 
We will show in the next section that our model is equivalent to the model of Ref. H for a given range of parameters. 



Fig. |9| represents the integration in the case (i > ^ (the importance of the value of (3 will be highlighted in Sec. fV) 



Ref. also reported an experiment in which the slider was stopped abruptly, but the applied stress was not 
released. Under these conditions, the medium does not recompactify towards the initial state but remains dilated 
in an intermediate state. This feature cannot be captured by our model, since the evolution of z does not explicitly 
depend on the applied stress but only on the horizontal velocity. In order to describe this effect, we modify Eq. (3) 
in order to explicitly include a stress dependence in the evolution of the dilatancy 

^ — 5—' y^) 



r] R 

where F^xt = k{Vt — x) is the applied force and A is a constant. The behavior of this model is similar to the 
simpler model introduced in Section II, but the zero velocity fixed point explicitly depend on the applied stress (i.e. 
z* = AFext)- Fig. Il^ shows the solution of the model compared whith the experiment of Ref. B. 



IV. LINEAR STABILITY 



The simple form of Eqs. (4-6) allows us to study analytically the linear stability of the system. We first concentrate 
on the inertial case and describe the main results about the dynamics of our problem (fixed point, critical curve). 
Next we discuss the origin of the instability and the connections with other models. Finally we investigate the nature 
of the bifurcation. 



A. Inertial case 



All the numerical results presented above have been obtained including inertial effects. The system of Eqs. (4-6) 
has a simple fixed point 

= '"'RkTWk' = ^' = "™]?TW 

We see that Zc tends to Zm when V tends to infinity, in agreement with experimental result. The critical line can also 
be computed explicitely in the framework of linear stability analysis. We skip the details of the calculations and just 
give the result 

^ ZraV^rjE? - Zmyq^Rfi + mR^u - mR^-qP + 2mR'^riViy — 2niRrfV (3 + mRTfV'^v - V'^rmfP 
" ,^7j^R^{R + f]V) ■ 
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ITll an d Fig. show the phase diagram in the k,V plane for the parameters used previously (in Sec. Ill A and 
III B| ). For both dry and wet granular layers we recover the stick-slip regime at sufficiently small k and V. In the 



Fig. 
Sec. 

dry case the critical velocity is higher than in the wet case and we can also identify the inertial regime on the right 
hand side of the phase diagram (see Fig. [ll| ). 



B. Non inertial case 



If we are interested only in low velocity displacements, the dynamical bifurcation line can be easily computed 
neglecting the mass of the slider 

Also in this case the dynamics is unstable for k below the critical line, but there is no inertial regime. We have no 
experimental results to compare with this relation which links all the relevant parameter of the model. 

Due to the simplicity of the non inertial case, we can write our system in the traditional form of a Hopf bifurcation 
p9[ , and calculate the coefficient which determines the nature of the transition. Without the inertial term this 
coefficient simply reduces to zero and therefore we have no information about the nature (super or subcritical) of 
the transition without pushing the caluclation to higher orders or including inertia. However, the calculation is 
particularly complex so we only analyze the problem mmrerically (see Sec. IVD). 



C. Dynamical friction force 



The stick-slip instability is due to the dependence of the friction coefficient on the velocity. Here we compute 
the friction force corresponding to the fixed point and show that the sign of P/ R — I'/rj plays an important role to 
determine the presence of an instability. In the steady state the friction force is given by 

Fc^Fd + Zrn jj^ .A P + l^V). (11) 

K + rjV 

For sufficiently high V, Fc does not depend on V, in agreement with experiments, but for relatively small velocities 
Fc{V) depends on V. The first derivative of the force is 

dF^V) ^ Rtj (3 V 
dV R + 7jV^R t]'' ^ ' 

We can thus identify three cases: 

• if {[3/ R — v/rj) is positive then there is a fc verifying Eq. (10) and below this k the system is unstable (the derivative 
in V of Fc{V) is negative i.e Fi.{V) decreases with V ). 

• if {/3/R — v/rf) is negative the system is always stable (fc cannot be negative). 

• if {(3/R — v/r] — Q) then Fc{V) does not depend on V. In this case the system is stable and we can write the 
friction force as 

F{z,x)=Fs + vz, (13) 



with Fs = Fd + vzm/rj. The form given in Eq. ( |13| ) for the friction force, together with Eq. (3), implies a friction 
coefficient independent on V and a stable steady state for all the values of the parameters. In the limit 77 3> 1, and 
assuming Eq. (p^), [3 tends to and the dilatancy rate is given by 



'X ^ , (14) 



which reproduces the model of Ref. 
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D. Nature of the bifurcation line 



The calculation in the non inertial case does not allow us to know the exact nature of the transition. Thus we 
investigate this problem numerically: the system is perturbed near its fixed point in the vertical position with different 
displacements (0.4/im — 0.002/Ltm). Two final states can be obtained depending on the position in the phase diagram: 
the system can evolve to the steady state or be driven to the stick-slip cycle. In an intermediate zone depending on 
the strength of the perturbation, the system can recover both the fixed point or the stick-slip regime. 

We identify three regimes, the first corresponds to the stick-slip regime where, independently on the amplitude of 
the perturbation the system falls into a periodic cycle. In the second regime, associated with high driving velocity, 
the system evolves to the stable fixed point. In the third intermediate regime, the final state depends on the initial 
perturbation: if the perturbation is sufficiently large the system falls into a periodic regime, while if the perturbation 
is weak it evolves towards the fixed point. The transition between the two regimes is discontinuous (i.e subcritical). 
Fig. |l^ shows the amplitude of the oscillations as a function of the driving velocity. It would be interesting to check 
experimentally the hysteretic nature of the bifurcation line. 

The presence of an hysteretic transition line could be related to an underlying first-order phase transition in the layer 
density induced by the applied stress. Recently analyzing the results of photoelastic disks in a two dimensional 
shear cell, it has been argued that the density of a the granular packing would be the order parameter of a second 
order phase transition induced by shear. It would be interesting to relate the different experimental phase transitions 
through a suitable microscopic model. 

V. DISCUSSION AND OPEN PROBLEMS 

We have introduced a model to describe the friction force of a sheared granular medium, treating explicitly the 
dilatancy during the slip, in the framework of rate and state constitutive equations. This approach allows us to include 
in the description the effect of the movement of the grains and the dependence of the friction coefficient on the dynamics 
of the layer. The variables used are mean-field like, since they represent macroscopic quantities like the position or the 
velocity but they are sufficient to describe phenomenologically the system. We have integrated the model for two sets 
of parameters, in order to make quantitative predictions for two different experimental configurations corresponding 
to dry and wet granular media. 

The results are in good agreement with experiments. In particular, we recover the hysteretic dependence of the 
friction force on velocity and obtain a good fit to the experimental data recorded in dry granular media. The effect of 
the weight of the slider plate is included in the model and allows us to recover a stick-slip instability at low velocity. 
The physical origin of the instability is then directly related to the recompaction of the material under normal stress. 
The dynamical phase diagram is calculated analytically both in the inertial and non inertial cases and inertia is found 
to change only the high velocity part of this diagram. The equations used to model the dependence of the friction law 
on the external parameters include explicitly the effect of recompaction in the evolution of the vertical slider position. 

The use of constitutive equations to model the friction force on complex interfaces is the simplest way to obtain 
quantitative results on the dynamics of the system. This approach provides good results in various fields, from 
geophysics or to nanotribology. In order to include the dynamics (or thermodynamics in the case of lubrication) of 
the material in the description, we need detailed informations about the material used. Our knowledge of sheared 
granular media is very poor due to the particulate and disordered nature of such materials and it is difficult to 
characterize the internal stress and strain rate. A precise description of the friction force for granular systems should 
include some information about the stress repartition inside the sheared material. This is a difficult problem which 
even for the simple case of a static pile cannot be solved completely. In the dynamical regime, the velocity depends on 
the precise nature of the contacts and on the friction force induced by them. Statistical models are needed to obtain 
a more complete macroscopic description based on the microscopic grain dynamics. In this respect, the analogy with 
phase transitions could be extremely fruitful. 

Experiments on granular flow over a rough inclined plane display an interesting behavior [ pT|p^ , which is ruled 
by frictional properties. The dynamic stops abruptly when the drag force decreases and the system freezes with 
grains remaining in a static configuration. These phenomena can be related to the dependence of the friction force 
on the velocity of the grains: an increase of the friction force when the velocity of the layer decreases can produce an 
instability as in the system discussed here. It will be interesting to see if the methods discussed in this paper can be 
applied to this and other situations. 

We thank J. S. Rice, and S. Roux for useful discussions and encouragements. We are grateful to J-C. Geminard for 
providing us with the data of his experiments and for interesting remarks. S. Z. is supported by EC TMR Research 
Network under contract ERBFMRXCT960062. 
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FIG. 1. Friction coefficient as a function of the slider ve- 
locity. The circles are the experimental data from Ref. [3] the 
line is the result of the model with parameters: u = %QKgs~^ , 
r] — 0.1s, m = ll.ig, R = 8/^m, Zm ~ Sfim, (3 — 3-^, 
k = 135iVm-^ 
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FIG. 2. Numerical calculation for dry granular me- 
dia in the stick-slip regime. iv — 90Kgs~^, rj = 0.1s, 
m = lOg, R = 8fim, 2„ = 8fim, f3 = 3^^, V = llOAtms"^ 
k = 135Ai'm~^). (a) spring elongation, (b) horizontal dis- 
placement, (c) slider velocity, (d) vertical displacement . 
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FIG. 3. Slider velocity for different loading velocities dur- 
ing a slip event, {u = 90Kgs~^, r) = 0.1s, m = lOg, R — 8fim, 
Zm = Sum, P = 3^, k = 135A'^m~^). Close to the transi- 
tion point the maximum velocity is weakly dependent on the 
driving velocity. 
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FIG. 4. Evolution to steady sliding state; numerical cal- 
culation for dry granular media, {u = 90Kgs~^, r/ = 0.1s, 
m = lOg, R — 8fim, Zm = 8fJ,m, (3 = 3-^, V = 500/ims~^, 
(a) spring elongation, (b) vertical position. 
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FIG. 5. 

m — lOg, R = Sum, 
k = 135Nm~^). The amplitude of the cycle depends on the 
initial condition. Here xq — Xc, vq = Vc, zq = Zc — e, where e 
is given in the legend, {xc, Vc, Zc) is the fixed point. 
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FIG. 6. The system is sheared at different constant veloci- 
ties. We show the effect of the velocity change on the spring 
elongation and on dilatancy, see Ref. [15-17]. {u = 90Kgs~^, 



rj = 0.1s, m — lOg, 
y = 50 500/ims-\ k 



R = 8/im, 



Zm = 8^m, /3 = 3- 
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FIG. 7. Numerical calculation for wet granular media in 
the stick-slip regime, {v = IQQOOK gs~^, rf = 100s, m = IQg, 
R = lOOfim, Zm = lO/xm, /3 = 15^, V = 0.2/ims-\ 
(a) spring elongation, (b) plate velocity, 
(c) vertical position. 
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FIG. 9. Steady sliding regime, effect due to (3 > 
{u = 10000ir3s"\ ri = 100s, m = Wg, R = lOO/im, 
Zm = W^im, (3 = 15^, V = 2/ims~\ k = llOiVm"^). (a) 
spring elongation, (b) horizontal position, (c) plate velocity. 
With (3 — and r; S> 1 we recover the result of Ref.[4]. 
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FIG. 8. Steady sliding regime, for the wet case. (3 — 
[v = 10000ii'3s"\ r) = 100s, m = IO5, 7? = 100/xm, 
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FIG. 10. The plate is stopped in the steady sliding regime 
(at time t ~ 60s) keeping the applied stress constant. We 
compare the solution of the model with the experimental data 
of Ref.[4]. {u = lOOOOK gs-^, r) = 30s, m = 14.5g, R = 59^im, 
Zm = 6^m, (3 = 0.01^, V = 28.35Atms-\ k = 189.5Nm-^, 
A = Q.137^miV-^). 
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FIG. 11. Dynamical phase diagram. Parameters here have 
been chosen in order to reproduce the data of Ref.[3]. on 
dry granular media, the circles are the results of numeri- 
cal calculation, — 90Kgs~^, rj = 0.1s, m — 11. Sg, 
R — 8fim, Zm = 8fj.m, /3 — 3-^). We can identify three 
different zones. At low velocity and low k the system exhibits 
a stick-slip instability, while at high velocity we observe an 
inertial regime. We also observe an intermediate region char- 
acterized by steady sliding. 
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FIG. 13. Amplitude of the oscillation as a function of driv- 
ing velocity. We perturb the system near its fixed point in 
the vertical direction. Two amplitudes are used in order 
to identify the three regimes, the diagram represents the fi- 
nal state of the system as a function of the perturbation e. 
(v = QQKgs~^, r) — 0.1s, m = lOg, R — 8/im, Zm — 8/im, 

p = 3^, k = mmm-^). 
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FIG. 12. Dynamical phase diagram for parameters cor- 
responding to the wet case. At high V and k we have 
a steady-sliding phase which becomes unstable at fc = fc*. 
{v = IQmQKgs'^, ri = 100s, m = Wg, R = lOO^m, 
Zm = 10/im, /3 = 15^). The circles represent numerical 
calculation. 
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